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High -2 galaxy redshift surveys open up exciting possibilities for precision determinations of neu¬ 
trino masses and inflationary models. The high -2 surveys are more useful for cosmology than low -2 
ones owing to much weaker non-linearities in matter clustering, redshift-space distortion and galaxy 
bias, which allows us to use the galaxy power spectrum down to the smaller spatial scales that are 
inaccessible by low -2 surveys. We can then utilize the two-dimensional information of the linear 
power spectrum in angular and redshift space to measure the scale-dependent suppression of matter 
clustering due to neutrino free-streaming as well as the shape of the primordial power spectrum. To 
illustrate capabilities of high -2 surveys for constraining neutrino masses and the primordial power 
spectrum, we compare three future redshift surveys covering 300 square degrees at 0.5 < 2 < 2, 
2 < 2 < 4, and 3.5 < 2 < 6.5. We find that, combined with the cosmic microwave background data 
expected from the Planck satellite, these surveys allow precision determination of the total neutrino 
mass with the projected errors of ct(to„, tot) = 0.059, 0.043, and 0.025 eV, respectively, thus yielding 
a positive detection of the neutrino mass rather than an upper limit, as tr(m„,tot) is smaller than the 
lower limits to the neutrino masses implied from the neutrino oscillation experiments, by up to a 
factor of 4 for the highest redshift survey. The accuracies of constraining the tilt and running index 
of the primordial power spectrum, o(n s ) = (3.8, 3.7, 3.0) x 10 -3 and a(a s ) = (5.9, 5.7, 2.4) x 1CU 3 
at ko = 0.05 Mpc -1 , respectively, are smaller than the current uncertainties by more than an or¬ 
der of magnitude, which will allow us to discriminate between candidate inflationary models. In 
particular, the error on a s from the future highest redshift survey is not very far away from the 
prediction of a class of simple inflationary models driven by a massive scalar field with self-coupling, 
a s = —(0.8 — 1.2) x 10 -3 . 

PACS numbers: 95.55.Vj,98.65.Dx,98.80.Cq,98.70.Vc,98.80.Es 


I. INTRODUCTION 

We are living in the golden age of cosmology. Vari¬ 
ous data sets from precision measurements of tempera¬ 
ture and polarization anisotropy in the cosmic microwave 
background (CMB) radiation as well as those of matter 
density fluctuations in the large-scale structure of the 
universe mapped by galaxy redshift surveys, Lyman-a 
forests and weak gravitational lensing observations are 
in a spectacular agreement with the concordance ACDM 
model mm- These results assure that theory of cos¬ 
mological linear perturbations is basically correct, and 
can accurately describe the evolution of photons, neu¬ 
trinos, baryons, and collisionless dark matter particles 
[5, fi 0 ) for given initial perturbations generated during 
inflation |&, 31 ■ The predictions from linear perturbation 
theory can be compared with the precision cosmological 
measurements, in order to derive stringent constraints on 
the various basic cosmological parameters. Future obser¬ 
vations with better sensitivity and higher precision will 
continue to further improve our understanding of the uni¬ 
verse. 

Fluctuations in different cosmic fluids (dark matter, 
photons, baryons, and neutrinos) imprint characteristic 
features in their power spectra, owing to their interac¬ 
tion properties, thermal history, equation of state, and 
speed of sound. A remarkable example is the acoustic 
oscillation in the photon-baryon fluid that was generated 
before the decoupling epoch of photons, 2 ~ 1088, which 
has been observed in the power spectrum of CMB tem¬ 


perature anisotropy [lO.], temperature-polarization cross 
correlation mi. and distribution of galaxies il- 

Yet, the latest observations have shown convincingly 
that we still do not understand much of the universe. The 
standard model of cosmology tells us that the universe 
has been dominated by four components. In chronolog¬ 
ical order the four components are: early dark energy 
(also known as “inflaton” fields), radiation, dark mat¬ 
ter, and late-time dark energy. The striking fact is that 
we do not understand the precise nature of three (dark 
matter, and early and late-time dark energy) out of the 
four components; thus, understanding the nature of these 
three dark components has been and will continue to be 
one of the most important topics in cosmology in next 
decades. Of which, one might be hopeful that the next 
generation particle accelerators such as the Large Hadron 
Collider (coming on-line in 2007) would find some hints 
for the nature of dark matter particles. On the other 
hand, the nature of late-time dark energy, which was dis¬ 
covered by measurements of luminosity distance out to 
distant Type la supernovae Bin, is a complete mys¬ 
tery, and many people have been trying to find a way to 
constrain properties of dark energy (see, e.g., [l6| for a 
review). 

How about the early dark energy, inflaton fields, which 
caused the expansion of the universe to accelerate in the 
very early universe? We know little about the nature 
of inflaton, just like we know little about the nature of 
late-time dark energy. The required property of infla¬ 
ton fields is basically the same as that of the late-time 
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dark energy component: both must have a large negative 
pressure which is less than —1/3 of their energy density. 
To proceed further, however, one needs more information 
from observations. Different inflation models make spe¬ 
cific predictions for the shape of the power spectrum 3] 
(see also Appendix B) as well as for other statistical prop¬ 
erties 0 of primordial perturbations. Therefore, one of 
the most promising ways to constrain the physics of in¬ 
flation, hence the nature of early dark energy in the uni¬ 
verse, is to determine the shape of the primordial power 
spectrum accurately from observations. For example, the 
CMB data from the Wilkinson Microwave Anisotropy 
Probe 0, combined with the large-scale structure data 
from the Two-Degree Field Galaxy Redshift Survey TSj, 
have already ruled out one of the popular inflationary 
models driven by a self-interacting massless scalar field 
M- Understanding the physics of inflation better will 
likely provide an important implication for late-time dark 
energy. 

“Radiation” in the universe at around the matter- 
radiation equality mainly consists of photons and neu¬ 
trinos; however, neutrinos actually stop being radiation 
when their mean energy per particle roughly equals the 
temperature of the universe. The physics of neutrinos 
has been revolutionized over the last decade by solar, 
atmospheric, reactor, and accelerator neutrino experi¬ 
ments having provided strong evidence for finite neutrino 
masses via mixing between different neutrino flavors, the 
so-called neutrino oscillations namsunm. These 
experiments are, however, only sensitive to mass square 
differences between neutrino mass eigenstates, implying 
Amjj ~ 7 x 10 -5 eV 2 and Am| 2 — 3 x 10 -3 eV 2 ; thus, 
the most fundamental quantity of neutrinos, the abso¬ 
lute mass, has not been determined yet. Cosmological 
neutrinos that are the relic of the cosmic thermal his¬ 
tory have distinct influences on the structure formation. 
Their large energy density, comparable to the energy den¬ 
sity of photons before the matter-radiation equality, de¬ 
termines the expansion history of the universe. Even 
after the matter-radiation equality, neutrinos having be¬ 
come non-relativistic affect the structure formation by 
suppressing the growth of matter density fluctuations at 
small spatial scales owing to their large velocity disper¬ 
sion [23, [Ui H3, EH ;29j, 30] (see Sec. II and Appendix A 
for more details). Therefore, the galaxy redshift surveys, 
combined with the CMB data, provide a powerful, albeit 
indirect, means to constraining the neutrino properties 
P EH likiL l34l Ip . This approach also complements 
the theoretical and direct experimental efforts for under¬ 
standing the neutrino physics. In fact, the cosmological 
constraints have placed the most stringent upper bound 
on the total neutrino mass, m^tot ~ 0.6 eV (2cr) _36$ . 
stronger than the direct experiment limit 2 eV |33 ■ In 
addition, the result obtained from the Liquid Scintillator 
Neutrino Detector (LSND) experiment, which implies 
to D e oscillations with Am 2 ^,0.2 eV 2 |38[ in an apparent 
contradiction with the other neutrino oscillation experi¬ 
ments mentioned above, potentially suggests the need for 


new physics: the cosmological observations will provide 
independent tests of this hypothesis. 

In this paper we shall study the capability of future 
galaxy surveys at high redshifts, combined with the CMB 
data, for constraining (1) the neutrino properties, more 
specifically the total neutrino mass, tot > and the num¬ 
ber of non-relativistic neutrino species, A" r , and (2) the 
shape of the primordial power spectrum that is parame¬ 
terized in terms of the spectral tilt, n s , and the running 
index, a s , motivated by inflationary predictions (see Ap¬ 
pendix B). For the former, we shall pay particular at¬ 
tention to our ability to simultaneously constrain m v ^ot 
and iV" r , as they will provide important clues to resolv¬ 
ing the absolute mass scale as well as the neutrino mass 
hierarchy. The accuracy of determining the neutrino pa¬ 
rameters and the power spectrum shape parameters will 
be derived using the Fisher information matrix formal¬ 
ism, including marginalization over the other cosmologi¬ 
cal parameters as well as the galaxy bias. 

Our analysis differs from the previous work on the 
neutrino parameters in that we fully take into account 
the two-dimensional nature of the galaxy power spec¬ 
trum in the line-of-sight and transverse directions, while 
the previous work used only spherically averaged, one- 
dimensional power spectra. The geometrical distortion 
due to cosmology and the redshift space distortion due to 
the peculiar velocity field will cause anisotropic features 
in the galaxy power spectrum. These features help to 
lift degeneracies between cosmological parameters, sub¬ 
stantially reducing the uncertainties in the parameter de¬ 
terminations. This is especially true when variations in 
parameters of interest cause modifications in the power 
spectrum shape, which is indeed the case for the neutrino 
parameters, tilt and running index. The usefulness of 
the two-dimensional power spectrum, especially for high- 
redshift galaxy surveys, has been carefully investigated 
in the context of the prospected constraints on late-time 
dark energy properties EH EH EH EH EH EH EHl • 

We shall show the parameter forecasts for future wide- 
held galaxy surveys that are already being planned or 
seriously under consideration: the Fiber Multiple Object 
Spectrograph (FMOS) on Subaru telescope EH; sig¬ 
nificantly expanded version, WFMOS [43, the Hobby- 
Ebery Telescope Dark Energy experiment (HETDEX) 
f48j ]. and the Cosmic Inflation Probe (CIP) mission EH- 
To model these surveys, we consider three hypothetical 
galaxy surveys which probe the universe over different 
ranges of redshift, (1) 0.5 < z < 2, (2) 2 < z < 4 and 
(3) 3.5 < 2 < 6.5. We hx the sky coverage of each sur¬ 
vey at Q s = 300 deg 2 in order to make a fair compari¬ 
son between different survey designs. As we shall show 
below, high-redshift surveys are extremely powerful for 
precision cosmology because they allow us to probe the 
linear power spectrum down to smaller length scales than 
surveys at low redshifts, protecting the cosmological in¬ 
formation against systematics due to non-linear pertur¬ 
bations. 

We shall also study how the parameter uncertainties 
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are affected by changes in the number density of sam¬ 
pled galaxies and the survey volume. The results would 
give us a good guidance to defining the optimal survey 
design to achieve the desired accuracies in parameter de¬ 
terminations. 

The structure of this paper is as follows. In Sec. un we 
review the physical pictures as to how the non-relativistic 
(massive) neutrinos lead to scale-dependent modifica¬ 
tions in the growth of mass clustering relative to the 
pure CDM model. Sec. HTT1 defines the parameterization 
of the primordial power spectrum motivated by inflation¬ 
ary predictions. In Sec. E3 we describe a methodology 
to model the galaxy power spectrum observable from a 
redshift survey that includes the two-dimensional nature 
in the line-of-sight and transverse directions. We then 
present the Fisher information matrix formalism that is 
used to estimate the projected uncertainties in the cos¬ 
mological parameter determination from statistical errors 
on the galaxy power spectrum measurement for a given 
survey. After survey parameters are defined in Sec. El 
we show the parameter forecasts in Sec. EH Finally, we 
present conclusions and some discussions in Sec. m We 
review the basic properties of cosmological neutrinos in 
Appendix A, the basic predictions from inflationary mod¬ 
els for the shape of the primordial power spectrum in Ap¬ 
pendix B, and the relation between the primordial power 
spectrum and the observed power spectrum of matter 
density fluctuations in Appendix C. 

In the following, we assume an adiabatic, cold dark 
matter (CDM) dominated cosmological model with flat 
geometry, which is supported by the WMAP results 
[Tj, 1 3(il l, and employ the the notation used in 0, M- 
the present-day density of CDM, baryons, and non- 
relativistic neutrinos, in units of the critical density, are 
denoted as H c , fib, and f l u , respectively. The total mat¬ 
ter density is then fl m = fl c + fib + fly, and /„ is the 
ratio of the massive neutrino density contribution to fl m : 
/y = fly/fl m . 


II. NEUTRINO EFFECT ON STRUCTURE 
FORMATION 


Structure formation is modified by non-relativistic 
neutrinos on scales below the Hubble horizon size 
when the neutrinos became non-relativistic, k nr = 
0.0145(my/leV) 1,/2 flm 2 liMpc -1 (see Eq. IA8I 1. In par¬ 
ticular, the characteristic scale imprinted onto the galaxy 
power spectrum at a given redshift z is the neutrino free- 
streaming scale, which is defined by Eq. E33J: 


k ls (z) = 0.113 Mpc 1 (j-^) 


/ Wy \ 

fn m h 2 

5 ^ 

VI eVJ 

V 0.14 

1 + zJ 


1/2 


( 2 ) 


Therefore, non-relativistic neutrinos with lighter masses 
suppress the growth of structure formation on larger spa¬ 
tial scales at a given redshift, and the free-streaming 
length becomes shorter at a lower redshift as neutrino 
velocity decreases with redshift. The most important 
property of the free-streaming scale is that it depends 
on the mass of each species, m„, rather than the total 
mass, N™m u ] thus, measurements of k{ s allow us to dis¬ 
tinguish different neutrino mass hierarchy models. For¬ 
tunately, fcf s appears on the scales that are accessible by 
galaxy surveys: kf s = 0.096 — 0.179 Mpc -1 at z = 6 — 1 
for my = 1 eV. 

On the spatial scales larger than the free-streaming 
length, k < k{ s , neutrinos can cluster and fall into gravi¬ 
tational potential well together with CDM and baryonic 
matter. In this case, perturbations in all matter com¬ 
ponents (CDM, baryon and neutrinos, denoted as ‘cb^’ 
hereafter) grow at the same rate given by 


D c bv{k,z) oc D(z) k<^kf s (z), (3) 


where D(z) is the usual linear growth factor (see, e.g., 
Eq. (4) in [Kj). On the other hand, on the scales smaller 
than the free-streaming length, k > kf s , perturbations in 
non-relativistic neutrinos are absent due to the large ve¬ 
locity dispersion. In this case, the gravitational potential 
well is supported only by CDM and baryonic matter, and 
the growth of matter perturbations is slowed down rela¬ 
tive to that on the larger scales. As a result, the matter 
power spectrum for k > kf s is suppressed relative to that 
for k < kfs . In this limit the total matter perturbations 
grow at the slower rate given by 


Throughout this paper we assume the standard ther¬ 
mal history in the early universe: there are three neutrino 
species with temperature equal to (4/11) 1 / 3 of the photon 
temperature. We then assume that 0 < lV“ r < 3 species 
are massive and could become non-relativistic by the 
present epoch, and those non-relativistic neutrinos have 
equal masses, m v . As we show in Appendix^ the den¬ 
sity parameter of the non-relativistic neutrinos is given 
by fly/i 2 = N™m v /(94.1 eV), where we have assumed 
2.725 K for the CMB temperature today [5(j, and h is the 
Hubble parameter defined as Hq = 100 h km s -1 Mpc -1 . 
The neutrino mass fraction is thus given by 


u 


Dy 

9m 


0.05 


\ 0.658 eV ) 



(1) 


D c \>v(k,z) oc (1 — /y)[D(^)] 1 p fc>fc fs (z), (4) 

where p = (5-^25 - 24/y)/4 |H. In [UIH an accurate 
fitting function for the scale-dependent growth rate was 
derived by matching these two asymptotic solutions. We 
shall use the fitting function throughout this paper. 

Figure E shows suppression in the growth rate of total 
matter perturbations at k = 0.01, 0.1, and 1 IiMpc -1 
due to the neutrino free-streaming. The suppression 
becomes more significant at lower redshifts for a given 
wavenumber, or for higher frequency perturbations at a 
given redshift, because neutrino can grow together with 
CDM and baryonic matter after the spatial scale of a 
given perturbation has become larger than the neutrino 
free-streaming scale that varies with redshift as given by 
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FIG. 1: Suppression in the growth rate of total matter per¬ 
turbations (CDM, baryons and non-relativistic neutrinos), 
D c b v (a), due to neutrino free-streaming. (a = (1 is the 

scale factor.) Upper panel. D c bv(a) / D v =o{a) for the neutrino 
mass fraction of f v = = 0.05. The number of non- 

relativistic neutrino species is varied from Nf r = 1, 2, and 3 
(from thick to thin lines), respectively. The solid, dashed, and 
dotted lines represent k = 0.01, 0.1, and 1 /iMpc -1 , respec¬ 
tively. Lower panel D c bv(a)/D„ = o(a) for a smaller neutrino 
mass fraction, /„ = 0.01. Note that the total mass of non- 
relativistic neutrinos is fixed to m v , tot = -N” r m„ = 0.66 eV 
and 0.13 eV in the upper and lower panels, respectively. 


Eq. (J2J). It is thus expected that a galaxy survey with 
different redshift slices can be used to efficiently extract 
the neutrino parameters, iV" r and m u . 

The upper and middle panels of Figure|21illustrate how 
free-streaming of non-relativistic neutrinos suppresses 
the amplitude of linear matter power spectrum, P(k), 
at z = 4. Note that we have normalized the primordial 
power spectrum such that all the power spectra match at 
k — * 0 (see § IIIIll . To illuminate the dependence of P(k) 
on to„, we fix the total mass of non-relativistic neutri¬ 
nos, N™m v , by /„ = 0.05 and 0.01 in the upper and 
middle panels, respectively, and vary the number of non- 
relativistic neutrino species as Nff = 1, 2 and 3. The 
suppression of power is clearly seen as one goes from 
k < k{ s (z ) to k > k{ s (z) (see Eq. [2] for the value of 
kf s ). The way the power is suppressed may be easily un¬ 
derstood by the dependence of k[ s (z) on to„; for example, 



FIG. 2: Upper panel. A fractional suppression of power in the 
linear power spectrum at z = 4 due to free-streaming of non- 
relativistic neutrinos. We fix the total mass of non-relativistic 
neutrinos by /„ = = 0.05, and vary the number of 

non-relativistic neutrino species (which have equal masses, 
m„) as JV° r = 1 (solid), 2 (dashed), and 3 (dot-dashed). 
The mass of individual neutrino species therefore varies as 
= 0.66, 0.33, and 0.22 eV, respectively (see Eq. 0). The 
shaded regions represent the 1-cr measurement errors on P(k) 
in each fc-bin, expected from a galaxy redshift survey observ¬ 
ing galaxies at 3.5 < z < 4.5 (see Table Q] for definition of the 
survey). Note that the errors are for the spherically averaged 
power spectrum over the shell of k in each bin. Different Nf r 
could be discriminated in this case. Middle panel. Same as 
in the upper panel, but for a smaller neutrino mass fraction, 
f v = 0.01. While it is not possible to discriminate between 
different Nf 1 , the overall suppression on small scales is clearly 
seen. Lower panel. Dependences of the shape of P{k) on the 
other cosmological parameters. 


P(k) at smaller k is more suppressed for a smaller m v , as 
lighter neutrinos have longer free-streaming lengths. On 
very small scales, k kf s (z) (k £ 1 and 0.1 Mpc -1 for 
fv = 0.05 and 0.01, respectively), however, the amount 
of suppression becomes nearly independent of k, and de¬ 
pends only on f v (or the total neutrino mass, N™m v ) 
as 


A P 

— 


2 /„ 


3 ln(H z=4 ) 
+ 5 


8 f v 


(5) 


We therefore conclude that one can extract f v and lV" r 
separately from the shape of P{k), if the suppression 
“pattern” in different regimes of k is accurately measured 
from observations. 
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Are observations good enough? The shaded boxes in 
the upper and middle panels in Figure 0 represent the 1- 
a measurement errors on P(k) expected from one of the 
fiducial galaxy surveys outlined in Sec. El We find that 
P(k) will be measured with ~ 1% accuracy in each k bin. 
If other cosmological parameters were perfectly known, 
the total mass of non-relativistic neutrinos as small as 
tot = A r " r TO y ^ 0.001 eV would be detected at more 
than 2 -a. This limit is much smaller than the lower 
mass limit implied from the neutrino oscillation exper¬ 
iments, 0.06 eV. This estimate is, of course, unrealistic 
because a combination of other cosmological parameters 
could mimic the N™ or f v dependence of P(k). The 
lower panel in Figure |2] illustrates how other cosmolog¬ 
ical parameters change the shape of P(k). In the fol¬ 
lowing, we shall extensively study how well future high- 
redshift galaxy surveys, combined with the cosmic mi¬ 
crowave background data, can determine the mass of non- 
relativistic neutrinos and discriminate between different 
A r “ r , fully taking into account degeneracies between cos¬ 
mological parameters. 

III. SHAPE OF PRIMORDIAL POWER 

SPECTRUM AND INFLATIONARY MODELS 

Inflation generally predicts that the primordial power 
spectrum of curvature perturbations is nearly scale- 
invariant. Different inflationary models make specific 
predictions for deviations of the primordial spectrum 
from a scale-invariant spectrum, and the deviation is of¬ 
ten parameterized by the “tilt”, n s , and the “running 
index”, a s , of the primordial power spectrum. As the pri¬ 
mordial power spectrum is nearly scale-invariant, \n s — 1| 
and |a s | are predicted to be much less than unity. 

This, however, does not mean that the observed mat¬ 
ter power spectrum is also nearly scale-invariant. In Ap¬ 
pendix [U| we derive the power spectrum of total matter 
perturbations that is normalized by the primordial cur¬ 
vature perturbation (see Eq. m) 


epoch. We use the fitting function found in [m} [3| for 
T(k). Note that the transfer function and the growth 
rate are normalized such that T(k ) —> 1 and D c b v /a —► 1 
as k —> 0 during the matter era. 

In Appendix m we describe generic predictions on n s 
and a s from inflationary models. For example, inflation 
driven by a massive, self-interacting scalar field predicts 
n s = 0.94 — 0.96 and a s = (0.8 — 1.2) x 10~ 3 for the num¬ 
ber of e-foldings of expansion factor before the end of 
inflation of 50. This example shows that precision deter¬ 
mination of n s and a s allows us to discriminate between 
candidate inflationary models (see Q for more details). 


IV. MODELING GALAXY POWER SPECTRUM 

A. Geometrical and Redshift-Space Distortion 

Suppose now that we have a redshift survey of galax¬ 
ies at some redshift. Galaxies are biased tracers of the 
underlying gravitational field, and the galaxy power spec¬ 
trum measures how clustering strength of galaxies varies 
as a function of 3-dimensional wavenumbers, k (or the 
inverse of 3-dimensional length scales). 

We do not measure the length scale directly in real 
space; rather, we measure (1) angular positions of galax¬ 
ies on the sky, and (2) radial positions of galaxies in 
redshift space. To convert (1) and (2) to positions in 
3-dimensional space, however, one needs to assume a ref¬ 
erence cosmological model, which might be different from 
the true cosmology. An incorrect mapping of observed 
angular and redshift positions to 3-dimensional positions 
produces a distortion in the measured power spectrum, 
known as the “geometrical distortion” [5J, |5J, [56} . The 
geometrical distortion can be described as follows. The 
comoving size of an object at redshift z in radial, rn, and 
transverse, rj_, directions are computed from the exten¬ 
sion in redshift, A z, and the angular size, Ad, respec¬ 
tively, as 


k 3 P(k, z) r2 f 2 k 2 \ 2 

(sflpd 

/ i \ — l+n s +^a„ ln(k/ko) 

xDl bv {k,z)T\k) , (6) 

where ko = 0.05 Mpc -1 , = 2.95 x 10 _9 A, and A 
is the normalization parameter given by the WMAP 
collaboration pj. We adopt A = 0.871, which gives 
S-jz = 5.07 x 10 -5 . (In the notation of jgjjjj, |§4} Sjz = 6^.) 
The linear transfer function, T(fc), describes the evolu¬ 
tion of the matter power spectrum during radiation era 
and the interaction between photons and baryons be¬ 
fore the decoupling of photons. Note that T{k) depends 
only on non-inflationary parameters such as f t m h 2 and 
and is independent of n s and a s . Also, the effects 
of non-relativistic neutrinos are captured in D c ^(k, z)', 
thus, T(k) is independent of time after the decoupling 


A2 

r|1 " H{ z y 

r± = D a (z)A9, 


(7) 


where D A is the comoving angular diameter distance 
given in the spatial sector of the Friedmann-Robertson- 
Walker line element, dl 2 = a 2 (dx 2 + D\dfl) (x is the 
comoving radial distance). We assume a flat universe 
throughout this paper, in which case \ = D A - The co¬ 
moving angular distance out to a galaxy at redshift z 
is 


D a (z) 



( 8 ) 


where H(z) is the Hubble parameter given by 
H 2 (z)=H 2 [n m (l + z) 3 + H a ] . 


( 9 ) 
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Here fi m + HA = 1, and Ha = A/(3Hq) is the present-day 
density parameter of a cosmological constant, A. A tricky 
part is that H(z) and D A {z) in Eq. 0 depend on cosmo¬ 
logical models. It is therefore necessary to assume some 
fiducial cosmological model to compute the conversion 
factors. In the following, quantities in the fiducial cos¬ 
mological model are distinguished by the subscript ‘fid’. 
Then, the length scales in Fourier space in radial, fegdii, 
and transverse, /cgdx, directions are estimated from the 
inverse of rfid || andrsax- These fiducial wavenumbers are 
related to the true wavenumbers by 


k± 


h 


Da (■£ )fid 

Da(z) 


fcfidX, 


H(z) 
H(z) fi d 


^fid|| • 


( 10 ) 


Therefore, any difference between the fiducial cosmolog¬ 
ical model and the true model would cause anisotropic 
distortions in the estimated power spectrum in (fcfidXj 
&fid||) space. 

In addition, shifts in 2 due to peculiar velocities of 
galaxies distort the shape of the power spectrum along 
the line-of-sight direction, which is known as the “redshift 
space distortion” j5jj . From azimuthal symmetry around 
the line-of-sight direction, which is valid when a distant- 
observer approximation holds, the linear power spectrum 
estimated in redshift space, P s {kf l <jj_, fegd||)> is modeled in 
as 


.Ps (fcfidX, &fid||) 


d a{z)I a H{z) 
D a (z) 2 H(z) fi d 

xb\P{k , z), 


1 + (3(k, z) 



( 11 ) 


where k = (k^_ + k p) 1 / 2 and 


P(k, z) 


1 din D c bis(k, z) 
b\ dln(l + 2) 


( 12 ) 


B. Comments on Baryonic Oscillations 

In this paper, we employ the linear transfer function 
with baryonic oscillations smoothed out (but includes 
non-relativistic neutrinos) M HI- As extensively in¬ 
vestigated in {3 |3 S3 , the baryonic oscillations can be 
used as a standard ruler, thereby allowing one to precisely 
constrain H(z) and Da{z) separately through the geo¬ 
metrical distortion effects (especially for a high-redshift 
survey). Therefore, our ignoring the baryonic oscillations 
might underestimate the true capability of redshift sur¬ 
veys for constraining cosmological parameters. 

We have found that the constraints on n s and a s from 
galaxy surveys improve by a factor of 2-3 when baryonic 
oscillations are included. This is because the baryonic os¬ 
cillations basically fix the values of H m , fl m h 2 and Hbh 2 , 
lifting parameter degeneracies between H m /i 2 , Hbh 2 , n s , 
and a s . However, we suspect that this is a rather opti¬ 
mistic forecast, as we are assuming a flat universe dom¬ 
inated by a cosmological constant. This might be a too 
strong prior, and relaxing our assumptions about geom¬ 
etry of the universe or the properties of dark energy will 
likely result in different forecasts for n s and a s . In this 
paper we try to separate the issues of non-flat universe 
and/or equation of state of dark energy from the physics 
of neutrinos and inflation. We do not include the bary¬ 
onic oscillations in our analysis, in order to avoid too 
optimistic conclusions about the constraints on the neu¬ 
trino parameters, n s , and a s . 

Eventually, the full analysis including non-flat uni¬ 
verse, arbitrary dark energy equation of state and its 
time dependence, non-relativistic neutrinos, n s , and a s , 
using all the information we have at hand including the 
baryonic oscillations, will be necessary. We leave it for 
a future publication (Takada and Komatsu, in prepara¬ 
tion). 


C. Parameter Forecast: Fisher Matrix Analysis 


is a function characterizing the linear redshift space 
distortion, and b\ is a scale-independent, linear bias 
parameter. Note that P(k,z) depends on both red¬ 
shift and wavenumber via the linear growth rate. In 
the infall regime, k <C kf s (z), we have bi(3(k,z) ~ 
—dlnD(z)/dln(l+z), while in the free-streaming regime, 
k fcf s (z), we have b\/3(k, z) « — (1— p)d\nD(z)/d\n(l+ 
z), where p is defined below Eq. Q. 

One might think that the geometrical and redshift- 
space distortion effects are somewhat degenerate in the 
measured power spectrum. This would be true only 
if the power spectrum was a simple power law. For¬ 
tunately, characteristic, non-power-law features in P(k) 
such as the broad peak from the matter-radiation equal¬ 
ity, scale-dependent suppression of power due to baryons 
and non-relativistic neutrinos, the tilt and running of 
the primordial power spectrum, the baryonic acoustic os¬ 
cillations, etc., hel p b reak degeneracies quite efficiently 

[3 El El El El El El El El 


In order to investigate how well one can constrain the 
cosmological parameters for a given redshift survey de¬ 
sign, one needs to specify measurement uncertainties of 
the galaxy power spectrum. When non-linearity is weak, 
it is reasonable to assume that observed density perturba¬ 
tions obey Gaussian statistics. In this case, there are two 
sources of statistical errors on a power spectrum measure¬ 
ment: the sampling variance (due to the limited number 
of independent wavenumbers sampled from a finite sur¬ 
vey volume) and the shot noise (due to the imperfect 
sampling of fluctuations by the finite number of galax¬ 
ies). To be more specific, the statistical error is given in 
Hby 


[A P a (*j)l 

2 2 

P s (ki) 

" N k [ 


1 + 


1 


n g P s {ki)_ 


(13) 


where n g is the mean number density of galaxies and N k 
is the number of independent k^d modes within a given 
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bin at /cfid = fc,: 


N k = 2nk 2 AkAp (^-^573 J • (14) 

Here 2ir/V s ' is the size of the fundamental cell in k 
space, 14 is the comoving survey volume, and p is the co¬ 
sine of the angle between fcfid and the line-of-sight. Note 
that we have assumed that the galaxy selection function 
is uniform over the redshift slice we consider and ignored 
any boundary effects of survey geometry for simplicity. 

The first term in Eq. CEl represents sampling vari¬ 
ance. Errors become independent of the number den¬ 
sity of galaxies when sampling variance dominates (i.e., 
P s n g over the range of k considered), and thus the 
only way to reduce the errors is to survey a larger vol¬ 
ume. On the other hand, the second term represents shot 
noise, which comes from discreteness of galaxy samples. 
When shot noise dominates (P B <C n g ), the most effec¬ 
tive way to reduce noise is to increase the number density 
of galaxies by increasing exposure time per field. Note 
that for a fixed n g the relative importance of shot noise 
contribution can be suppressed by using galaxies with 
larger bias parameters, b x , as P s oc b\. In Sec. |V| we shall 
discuss more about the survey design that is required to 
attain the desired parameter accuracy. 

We use the Fisher information matrix formalism to 
convert the errors on P s {k) into error estimates of model 
parameters l.fil . The Fisher matrix is computed from 


Pa/3 


A /■*“>* 

8 tH J_ x J ktnin dp a dp 0 


n g P s (k,p ) 
n g P s (k,p) +1 


(15) 


where p a expresses a set of parameters. One may evalu¬ 
ate some derivative terms analytically: 


d In P s (k,p) 
dS-ji 

dlnP s (k, p) 
dn s 

dlnP s (k,p) 

da s 


2 

5-r 1 

(16) 

1 k 
ln I“’ 
ko 

(17) 

Ifi„A) 2 . 

2 V k 0 J 

(18) 


The la error on p a marginalized over the other parame¬ 
ters is given by a 2 (p a ) = {F~ 1 ) aa , where F ~ 4 is the in¬ 
verse of the Fisher matrix. It is sometimes useful to con¬ 
sider projected constraints in a two-parameter subspace 
to see how two parameters are correlated. We follow the 
method described around Eq. (37) in |53| for doing this. 
Another quantity to describe degeneracies between given 
two parameters, p g and p„, is the correlation coefficient 
defined as 


r(p a ,pp) 


(F-'U 

x/(F-i) aa (F-i)w 


(19) 


If |r| = 1, the parameters are totally degenerate, while 
r = 0 means they are uncorrelated. 

To calculate F a p using Eq. G3» we need to specify 
fcmin and fc max for a given galaxy survey. We use the 
upper limit, fc max , to exclude information in the non¬ 
linear regime, where the linear theory prediction of den¬ 
sity fluctuations, Eq. 03- becomes invalid. Following 
[ 33 , we adopt a conservative estimate for fc max by im¬ 
posing the condition a mass (R, z) = 0.5, where a mass (R , z) 
is the r.m.s. mass fluctuation in a sphere of radius 
R = 7 t /(2 fc max ) at a given redshift z. All the Fourier 
modes below fc max are considered as in the linear regime. 
This idea is partly supported by the simulation-based 
work in the literature |60l . IfiTL l62il . while a more care¬ 
ful and quantitative study is needed to understand the 
impact of non-linearities on cosmological parameter esti¬ 
mates as well as to study how to protect the cosmological 
information against the systematics. Table [I] lists /c max 
for each redshift slice of galaxy surveys we shall consider. 
In addition, we shall show how the results will change 
with varying fc max . As for the minimum wavenumber, 
we use k min = 10 -4 Mpc -1 , which gives well-converged 
results for all the cases we consider. 


D. Model Parameters 

The parameter forecasts derived from the Fisher infor¬ 
mation formalism depend on the fiducial model and are 
also sensitive to the choice of free parameters. We include 
a fairly broad range of the CDM dominated cosmology: 
the density parameters are Q m (= 0.27), fl m h 2 (= 0.14), 
and flbh 2 (= 0.024) (note that we assume a flat universe); 
the primordial power spectrum shape parameters are the 
spectral tilt, n s (= 1 ), the running index, a s (= 0 ), and 
the normalization of primordial curvature perturbation, 
Sn(= 5.07 x 10 -5 ) (the numbers in the parentheses de¬ 
note the values of the fiducial model). The linear bias 
parameters, b x , are calculated for each redshift slice as 
given in Table [T| the fiducial values of the neutrino pa¬ 
rameters, f v and ./V" 1 ', are allowed to vary in order to 
study how the constraints on /„ and N™ change with 
the assumed fiducial values. For a survey which consists 
of N s redshift slices, we have 8 + N s parameters in total. 

As we shall show later, a galaxy survey alone can¬ 
not determine all the cosmological parameters simulta¬ 
neously, but would leave some parameter combinations 
degenerated. This is especially true when non-relativistic 
neutrinos are added. Therefore, it is desirable to com¬ 
bine the galaxy survey constraints with the constraints 
from CMB temperature and polarization anisotropy, in 
order to lift parameter degeneracies. When computing 
the Fisher matrix of CMB, we employ 7 parameters: 6 
parameters (the parameters above minus the neutrino 
parameters and the bias parameters) plus the Thom¬ 
son scattering optical depth to the last scattering sur¬ 
face, r(= 0.16). Note that we ignore the effects of non- 
relativistic neutrinos on the CMB power spectra: their 













effects are small and do not add very much to the con¬ 
straints from the high-2 galaxy survey. We then add 
the CMB Fisher matrix to the galaxy Fisher matrix as 
F af 3 = + F™ B . We entirely ignore the contri¬ 

bution to the CMB from the primordial gravitational 
waves. We use the publicly-available CMBFAST code 
jj65| to compute the angular power spectrum of tempera¬ 
ture anisotropy, C) TT , .E-mode polarization, C EE , and 
their cross correlation, C ; TE . Specifically we consider 
the noise per pixel and the angular resolution of the 
Planck experiment that were assumed in ’66]. Note that 
we use the CMB information in the range of multipole 
10 <l< 2000 . 

V. GALAXY SURVEY PARAMETERS 

We define the parameters of our hypothetical galaxy 
surveys so that the parameters resemble the future sur¬ 
veys that are already being planned and seriously pur¬ 
sued. As shown in Eq. Cl, the statistical error of the 
galaxy power spectrum measurement is limited by the 
survey volume, U s , as well as the mean number density 
of galaxies, n g . There are two advantages for the high- 
redshift galaxy surveys. First, given a fixed solid angle, 
the comoving volume in which we can observe galaxies is 
larger at higher redshifts than in the local universe. Ac¬ 
cordingly, it would be relatively easy to obtain the well- 
behaved survey geometry, e.g., a cubic geometry that 
would be helpful to handle the systematics. Second, den¬ 
sity fluctuations at smaller spatial scales are still in the 
linear regime or only in the weakly non-linear regime at 
higher redshift, which gives us more leverages on mea¬ 
suring the shape of the linear power spectrum. 

Of course, we do not always win by going to higher 
redshifts. Detecting galaxies at higher redshifts is obvi¬ 
ously more observationally demanding, as deeper imag¬ 
ing capabilities and better sensitivity for spectrographs 
are required. To increase the survey efficiency, the use 
of Multi Object Spectrographs (MOS) or Integral Field 
Unit (IFU) spectrographs will be favorable. It is there¬ 
fore unavoidable to have a trade-off in the survey design 
between the number of spectroscopic targets and the sur¬ 
vey volume: for a fixed duration of the survey (or a fixed 
amount of budget), the total number of spectroscopic 
targets would be anti-correlated with the survey volume. 
As carefully discussed in a survey having n g P g A 3 
over the range of wavenumbers considered is close to an 
optimal design. 

To make the comparison between different survey de¬ 
signs easier, we shall fix the total sky coverage of the 
surveys to 

survey ~ 300 ’leg - , 

for all cases. We choose to work with three surveys ob¬ 
serving at three different redshift ranges: 

• Gl: 0.5 < 2 < 2 


• G2: 2 < 2 < 4 

• SG: 3.5 < 2 < 6.5 

where Gl and G2 stands for the “Ground-based galaxy 
survey” 1 and 2, respectively, while SG stands for the 
“Space-based Galaxy survey”. Table 1 lists detailed sur¬ 
vey parameters for each survey design. 


A. Gl: Ground-based Galaxy Survey at 0.5 < 2 < 2 

The first survey design, Gl, is limited to 0.5 < 2 < 2 
for the following reasons. One of the reasonable target 
galaxies from the ground would be giant ellipticals or 
star-forming galaxies because of their large luminosity. If 
spectroscopic observations in optical wavebands are avail¬ 
able, galaxies having either 3727 A[OII] emission lines or 
4000 A continuum break may be selected, in which case 
2 = 1.3 would be the highest redshift bin, as these spec¬ 
tral line features will move out of the optical wavebands 
otherwise. If spectroscopy in the near infrared band is 
available, such as that proposed by the FMOS instru¬ 
ment on Subaru telescope, one may select 6563 A Ha 
emission lines which usually have the highest equivalent 
width among the lines in a star-forming galaxy, in which 
case a higher redshift, 2 A 2, may be reached. Based on 
these considerations, we consider a survey of 0.5 < 2 < 2 
and subdivide the survey into 3 redshift bins centering at 
2 = 0.75, 1.25 and 1.75 with widths A2 = 0.5. While it 
is currently difficult to estimate the number density and 
bias parameters for these galaxies with any certainty be¬ 
cause we have a limited knowledge of how such galaxies 
formed within the CDM hierarchical clustering scenario, 
we follow the argument given in mmm and assume 
n g = 0.5 x 10 -3 h 3 Mpc -3 . We determine b\ so that it 
satisfies the condition as^ g = 1 at a given redshift, where 



with /3 m = —dlnD/d\n(l + 2) (i.e. we do not include 
the massive neutrino contribution when we estimate the 
fiducial b\). Note that there is no a priori reason to be¬ 
lieve that the r.m.s. fluctuation of the number density of 
galaxies within an 8 h ~ 1 Mpc sphere should be unity; this 
condition is rather motivated by observations, and it does 
not have to be true for arbitrary population of galaxies. 
Nevertheless, this approach seems to provide reasonable 
values for 61, and also it makes it easier to compare our 
results with the previous work that used the same recipe 
mm. The values of fc max in Table [I] are computed by 
o’mass(-R, 2) = 0.5, where cr mass (i?, 2) is the r.m.s. mass 
fluctuation in a sphere of radius corresponding to fc max , 
R = 7r/(2fc max ) (see Sec. II V All . 
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Survey 

-^center 

^max 

(/iMpc -1 ) 

^survey 

(deg 2 ) 

U 

{h~ 3 Gpc 3 ) 

n B 

(10 -3 /i 3 Mpc -3 ) 

Bias 

PgTlg 

(^max) 

Gl (0.5 < 2 < 2) 

0.75 

0.14 

300 

0.33 

0.5 

1.22 

4.83 


1.25 

0.19 

300 

0.53 

0.5 

1.47 

2.49 


1.75 

0.25 

300 

0.64 

0.5 

1.75 

1.38 

G2 (2 < 2 < 4) 

2.25 

0.32 

300 

0.68 

0.5 

2.03 

0.80 


2.75 

0.41 

300 

0.69 

0.5 

2.32 

0.46 


3.25 

0.52 

300 

0.67 

0.5 

2.62 

0.27 


3.75 

0.64 

300 

0.64 

0.5 

2.92 

0.16 

SG (3.5 < 2 < 6.5) 

4 

0.71 

300 

1.26 

5 

4 

2.19 


5 

1.01 

300 

1.13 

5 

5 

1.04 


6 

1.50 

300 

1.02 

5 

5.5 

0.35 


TABLE I: Galaxy survey specifications that we assume in this paper (see Sec. El for the details). We assume a fixed sky 
coverage (300 deg 2 ) for all the surveys, and V B and n g are the comoving survey volume and the comoving number density of 
sampled galaxies for each redshift slice, respectively, ^center denotes the center redshift of each redshift slice, and fc max is the 
maximum wavenumber below which information in the linear power spectrum can be extracted. (We do not use any information 
above fc max in the Fisher information matrix analysis.) “Bias” denotes the assumed linear bias parameters of sampled galaxies. 


B. G2: Ground-based Galaxy Survey at 2 < z < 4 

The second design, G2, probes higher redshifts than Gl 
by observing different tracers. The primary candidates 
from the ground in this redshift range would be Lyman- 
break galaxies or Lyman-a emitters, which are accessi¬ 
ble from a deep survey of 8-m class telescopes in optical 
wavebands. This type of survey has been proposed by the 
Hobby-Ebery Telescope Dark Energy experiment (HET- 
DEX) and Wide Field Multiple Object Spectrograph 
(WFMOS) collaborations (42$ • To make the comparison 
easier, we shall assume the same number density as for 
Gl, 0.5 x 10 -3 h 3 Mpc -3 , for G2. This number cor¬ 
responds to 4500 galaxies per square degrees, or 1.25 
galaxies per square arcminutes for the surface density. 
We subdivide the redshift range of G2 into 4 bins, 2.25, 
2.75, 3.25 and 3.75, and again determine the bias param¬ 
eters by imposing as i9 (z) = 1 at center redshifts of each 
bin. 


C. SG: Space-based Galaxy Survey at 4 < a < 6 

The third design, SG, is a space-based observation 
which targets galaxies at even higher redshifts, 4 ^ 6. 

The useful line features will be redshifted into infrared, 
which makes such high -2 galaxies accessible only from 
space. We determine the survey parameters on the ba¬ 
sis of the Cosmic Inflation Probe (CIP) mission [ 4 ^, one 
of the nine studies selected by NASA to investigate new 
ideas for future mission concepts within its Astronomical 
Search for Origins Program. The CIP is a slitless-grating 
survey in the near infrared, 2.5-5 /xm, which detects Ha 
emission lines in star-forming galaxies at these redshifts. 
Being up in space with low background, CIP can achieve 
a superb sensitivity in infrared. We assume the num¬ 


ber density of 5 x 10 -3 h 3 Mpc -3 4Sj, which is larger 
by one order magnitude than that by the ground based 
surveys. This number density may be partly justified by 
the fact that the Lyman break galaxies or the Lyman-a 
emitters show the similar number density at these red¬ 
shifts as implied from a deep imaging survey Ifi7ll . and 
most of such galaxies are very likely to exhibit an even 
stronger Ha emission line. For this survey, we assume 
the bias parameters, b\ = 4.5, 5 and 5.5 for redshift 
slices of 2 = 4, 5 and 6 with redshift width A 2 = 1, 
respectively. The bias parameters for this survey have 
been determined using a different method. We used the 
mass-weighted mean halo bias above a certain minimum 
mass, M m ; n . The minimum mass was found such that 
the number density of dark matter halos above Af m ; n 
should match the assumed number density of galaxies, 
JdM dn/dM = 5 x 10 -3 h 3 Mpc -3 . Therefore, we 
basically assumed that each dark matter halo above M m ; n 
hosts one Ha emitter on average. One may improve this 
model by using the Halo Occupation Distribution model, 
at the expense of increasing the number of free parame¬ 
ters. 

Note that we chose these survey designs not to say 
these are the optimal designs for doing cosmology with 
high -2 surveys, but rather to show how well these planned 
surveys can constrain the neutrino and inflationary pa¬ 
rameters. We are hoping that our results provide some 
useful information in designing high -2 galaxy surveys. 


VI. PARAMETER FORECAST: BASIC 
RESULTS 

Tables El and EH summarize the basic results of our 
forecasts for the cosmological parameters from the high- 
2 galaxy redshift surveys combined with the Planck data. 
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Survey 

f v {rn u , tot eV) 

Km = 
N™ 

1 A“ fid = 
AT 

2 A” r fid = 3 
AT 

n s 

a s 

111 

In 5n 

In I2 m /j 2 

In 12b h 2 

Planck alone 

- 

- 

- 

- 

0.0062 

0.0067 

0.035 

0.013 

0.028 

0.011 

G1 

0.0045(0.059) 

0.31 

0.64 

1.1 

0.0038 

0.0059 

0.0072 

0.0099 

0.0089 

0.0075 

G2 

0.0033(0.043) 

0.20 

0.49 

0.90 

0.0037 

0.0057 

0.0069 

0.0099 

0.0086 

0.0072 

SG 

0.0019(0.025) 

0.14 

0.40 

0.80 

0.0030 

0.0024 

0.0041 

0.0090 

0.0055 

0.0050 

All (G1+G2+SG) 

0.0018(0.024) 

0.091 

0.31 

0.60 

0.0026 

0.0023 

0.0030 

0.0089 

0.0043 

0.0048 


TABLE II: The projected 68% error on the cosmological parameters from Planck’s CMB data alone (the 1st row) and the 
high-z galaxy survey data combined with the Planck data (from the 2nd to 5th rows). The quoted error for a given parameter 
includes marginalization over the other parameter uncertainties. Note that the values with and without parenthesis in the 2nd 
column are the errors for /„ and m„,tot (eV), respectively. The fiducial values for the neutrino parameters are fv,Rd = 0.05 
(the neutrino mass fraction; Eq. 0) and A“ r fid = 3 (the number of non-relativistic neutrinos), while the fiducial values for the 
other parameters are given in Sec. EYE We also vary the fiducial value of from 1 to 2 and 3 when quoting the projected 
errors for A” r with fv,Rd = 0.05 being fixed, as indicated in the 3rd and 5th columns. 


Survey 

fv{mv, tot eV) 

A^d = 
A“ r 

1 A„,fid = 
A” r 

2 A^fid = 3 
A? r 

Us 

a s 

Om 

In 5k 

In Q.mh 2 

In 12b h 2 

G1 

0.0044(0.058) 

2.2 

7.1 

14 

0.0037 

0.0059 

0.0069 

0.0099 

0.0085 

0.0073 

G2 

0.0033(0.043) 

2.1 

7.1 

13 

0.0036 

0.0058 

0.0059 

0.0098 

0.0075 

0.0066 

SG 

0.0021(0.028) 

1.9 

6.4 

13 

0.0028 

0.0021 

0.0034 

0.0090 

0.0048 

0.0048 

All (G1+G2+SG) 

0.0019(0.025) 

1.1 

3.7 

7.4 

0.0021 

0.0016 

0.0017 

0.0087 

0.0030 

0.0045 


TABLE III: Same as in the previous table, but for the smaller fiducial neutrino mass fraction, fv,Rd = 0.01. 


Each column shows the projected 1 -ct error on a partic¬ 
ular parameter, marginalized over the other parameter 
uncertainties. The 1st row in Table EH shows the con¬ 
straints from the Planck data alone, while the other rows 
show the constraints from the Planck data combined with 
each of the high-z galaxy surveys outlined in the previ¬ 
ous section. The final row shows the constraints from 
all the data combined. The difference between these two 
Tables is the fiducial value for /„ = Q, ,/f2 m : Table ITTI 
uses f v = 0.05, whereas Table EH uses a lower value, 
f v = 0.01, as the fiducial value. In addition, in each Ta¬ 
ble the fiducial value for the number of non-relativistic 
neutrino species, Aj) r , is also varied from A" r = 1 to 2 
to 3. Therefore, in Table EED the fiducial mass of indi¬ 
vidual non-relativistic species changes from rriv = 0.66 
to 0.33 to 0.22 eV, whereas in Table Hill it changes from 
m v = 0.13 to 0.066 to 0.044 eV for N™ = 1, 2, and 3, 
respectively. It is also worth showing how the parameter 
errors are correlated with each other, and we give the 
parameter correlations in Table EE for the case of SG 
combined with Planck in Appendix IT)1 


A. Neutrino Parameters 

In this paper, we are particularly interested in the ca¬ 
pability of future high-z redshift surveys to constraining 
the neutrino parameters, f v and A" r , as well as the shape 
of the primordial power spectrum, the tilt ( n s ) and the 


running spectral index (a s ). First, we study the neutrino 
parameters. 

The upper panel of Figure 0 shows error ellipses in the 
subspace of ( f v , A" r ). Two “islands” show two different 
fiducial models: the left island is (fv,N™) = (0.01,1), 
while the right island is (0.05,3). We find that the er¬ 
rors on f v and N™ are only weakly degenerate with each 
other, implying that the constraints on the two parame¬ 
ters come from different regions of P(k) in fc-space, which 
can be seen more clearly from Figure E 


1. Total Neutrino Mass 

As we have shown in Sec. II, galaxy surveys constrain 
the total mass of non-relativistic neutrinos by measuring 
the overall suppression of power at small scales compared 
with the scales larger than the neutrino free-streaming 
length, A P(k)/P(k) ~ —8 f v [Eq. ©]. Tables UTI and 
EH and Figure E (the widths of the error ellipses show 
the accuracy of constraining the total neutrino mass, 
A" r ?7i„) show that the high-z galaxy surveys can pro¬ 
vide very tight constraints on the total neutrino mass, 
mv ,tot = A“ r m„, and the constraint improves steadily 
by going to higher redshifts for a given survey area. 
The projected error (assuming f v = 0.05) improves from 
o^m^tot) = 0.059 to 0.025 eV for G1 to SG. (The con¬ 
straint on m J/j tot is very similar for f v = 0.01.) This 
is because the suppression rate of the amplitude of the 













11 


Vn^ [eV] 

0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 


GO 

0 

’o 

0 3 

Oh 

GO 

O 

G 


5 2 
0 
a 


P5 


* 1 


' 0.5 


1 1 1 i 1 1 1 1 i 1 1 1 1 i 1 | 1 1 i | 1 1 1 i 1 1 1 1 i 

N.O. exp. lower limit 


mal’ 

"inv< 


erse' m 3 =0,m 1 =m 2 



L ' N e „ 


, nnlM tvA d=0° 5 'N, >lid =3 

■A/,fid 0-01,N,.fid 1 

fi =0.27, Q m h 2 =0. 14 

m m 

SG+Planck 
SG+Planck (a=0) 
All+Planck (a=0) 

u— i — i — | — i —i—^~i— 1 — i — i — i —— i — i — i — i — I—i — i — i — i 



0 0.01 0.02 0.03 0.04 0.05 0.06 

f i/ = V n m) 


FIG. 3: Upper panel: Projected 68 % error ellipses in the 
neutrino parameter, (/„-A“ r ) plane, expected from the high- 
z galaxy survey data combined with the Planck data (see 
Table [I] for the survey definition). The two fiducial mod¬ 
els for f„ and N™ are considered: the left contours assume 
(.A, fid, = (0-01,1), while the right contours assume 

(A,fid, IV“g d ) = (0.05, 3). The outer thin lines and the middle 
light-gray contours are the forecasts for SG (the space-based 
mission at 3.5 < 2 < 6.5) plus Planck, without and with 
a prior on the running spectral index, a s = 0, respectively. 
The innermost, dark gray contours show the forecasts when 
all the galaxy surveys (two ground-based surveys and SG) and 
Planck are combined. The vertical dashed and dotted lines 
show the lower limits on /„ implied from the neutrino oscilla¬ 
tion experiments assuming the normal and inverted mass hi¬ 
erarchy models, respectively. The dashed and dotted curves 
then show the effective number of non-relativistic neutrino 
species, N e g, for the two hierarchy models (see Eq. 1211 for 
the definition). Lower panel: The projected 68 % on JV“ r as a 
function of the fiducial value of /„. The thick solid, dashed, 
and dotted lines use the fiducial values of N™\ d = 1, 2, and 3, 
respectively. The dot-dashed curve shows the difference be¬ 
tween JVeff for the normal and inverted mass hierarchy models. 
The leftmost thin solid line shows the error expected from a 
hypothetical full-sky SG survey for = 1. 


linear power spectrum on small scales can be precisely 
measured by the galaxy survey when combined with the 
tight constraint on the amplitude of the spectrum on 
large scales from Planck (also see Tablelvland Table IWfl 
in Appendix EJ). The steady improvement at higher red- 
shifts is simply because surveys at higher redshifts can 
be used to probe smaller spatial scales (i.e., larger fc max ; 
see Table I). We find that it is crucial to increase /e max 
as much as possible in order to improve the constraints 
on the neutrino parameters. Figure@]shows how one can 



FIG. 4: The projected 68 % error on f v (upper panel) and 
AT r (lower panel) against the maximum wavenumber, fc max , 
assuming the information in the linear power spectrum at k < 
fcmax can be used in the Fisher matrix analysis. The arrows 
in the above £-axis indicate our nominal fc max used in the 
analysis at each redshift (see Tablc0. Note that /„,Rd = 0.05 
and AT“fl d = 3 are assumed. 


reduce the errors on /„ and N™ by increasing fc max . 

These results are extremely encouraging. If we a pri¬ 
ori assume three-flavor neutrinos in compatible with the 
neutrino oscillation experiments, then it is likely that one 
can determine the sum of neutrino masses using the high¬ 
er galaxy redshift surveys, as the current lower bounds 
implied from the neutrino oscillation experiments are 
m„.tot ~ 0.06 and 0.1 eV for the normal and inverse mass 
hierarchies, respectively. It should also be noted that a 
detection of the total mass in the range of m„ )to t <0.1 eV 
gives an indirect evidence for the normal mass hierarchy, 
thereby resolving the mass hierarchy problem. 

Our results may be compared with the previous work 
[3l| . where a(m„ ito t) ~ 0.3 eV (lcr) was obtained; there¬ 
fore, our errors are smaller than theirs by a factor of 5-10, 
even though the survey volume that we assumed is larger 
only by a factor of 1.5-3 than what they assumed. What 
drives the improvement? There are two reasons. The 
first reason is because we consider high -2 galaxy surveys, 
while jEl considered low -2 surveys, such as the Sloan 
Digital Sky Survey, which suffer from much stronger non¬ 
linearity. We are therefore using the information on the 
power spectrum down to larger wavenumbers (i.e., k max 
is larger). Figure Q] shows that all the galaxy surveys 
have essentially equal power of constraining the neutrino 
parameters, when the information up to the same /c max 
is used. However, one cannot do this for low -2 surveys 
because of strong non-linearity. As long as we restrict 
ourselves to the linear regime, a higher redshift survey is 
more powerful in terms of constraining the neutrino pa¬ 
rameters. Interestingly, the error on the neutrino param¬ 
eters appears to be saturated at fc max ~ 1 Mpc -1 ; thus, 
the space-based mission, SG, is already nearly optimal for 
constraining the neutrino parameters for the survey pa- 


) [eV] 
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rameters (especially the number density and bias param¬ 
eters of sampled galaxies that determine the shot noise 
contribution to limit the small-scale measurements). The 
second reason is because our parameter forecast uses the 
full 2D information in the redshift-space power spectrum 
(see Eq. JTTI ) that includes effects of the cosmological 
distortion and the redshift-space distortion due to pecu¬ 
liar velocity. These effects are very useful in breaking 
parameter degeneracies. Table El shows that the con¬ 
straint on tot would be significantly degraded if we 
did not include the distortion effects. In particular, ig¬ 
noring the information on the redshift-space distortion, 
which is consistent with the analysis of (U, leads to a 
similar-level constraint on m ^,tot as theirs. The inclusion 
of the redshift-space distortion helps break degeneracy 
between the power spectrum amplitude and the galaxy 
bias, which in turn helps determine the small-scale sup¬ 
pression due to the neutrino free-streaming. 

The projected error on the neutrino total mass might 
also depend on the fiducial value of fl m h 2 , as the ef¬ 
fect of neutrinos on the powe r sp ectrum depends on 
fv oc ??Vtot/(^mft 2 ) [Eq. GJ] HU]. For a given /„, a 
variation in fl m h 2 changes m v ^ot- Now that fI m /i 2 has 
been constrained accurately by WMAP, however, we find 
that our results are not very sensitive to the precise value 
of fl m h 2 . We have repeated our analysis for S2 m h 2 =0.1, 
the 2 <t- level lower bound from the WMAP results Q|, and 
found very similar results. 


2. Number of Non-relativistic Neutrino Species 


Galaxy surveys could also be used to determine the 
individual mass of non-relativistic neutrinos, m v . As 
we have shown in Sec. II, galaxy surveys can constrain 
by determining the free-streaming scale, k{ s (z) oc m„ 
[Eq. ©], from distortion of the shape of the galaxy power 
spectrum near kf s (z). 

Neutrino oscillation experiments have provided tight 
limits on the mass square differences between neutrino 
mass eigenstates as | m 2 — m 2 | ~ 7 x 10 -5 eV 2 and 
[to 2 — m 2 | ~ 3 x 10 -3 eV 2 where to,; denotes mass of 
the I-tli mass eigenstates. We model a family of possi¬ 
ble models by the largest neutrino mass, m„, motivated 
by the fact that structure formation is sensitive to most 
massive species. We may define the effective number of 
non-relativistic neutrino species as 


v eff 


rrii rrij 
m v m v ’ 


( 21 ) 


which continuously varies between 1 < TV"g < 3. The 
total neutrino mass is given by m V} tot = AP) \m v . We 
then consider two neutrino mass hierarchy models (we 
shall assume, by convention, that m 2 > mi): 


• Normal mass hierarchy: m v = m 3 

In this model, m 3 is assumed to be the largest mass, 
and mi < m 2 < m 3 are allowed. When m 2 , mi <C 


m 3 (the extreme case mi = 0 eV), TV"^ ~ 1. On 
the other hand, when mi and m 2 are comparable 
to the mass difference between m 2 and m 3 , we have 
N^g ~ 3 (i.e., three masses are nearly degenerate). 
Hence, the “normal mass hierarchy” model allows 
Nfg to vary in the full parameter space, 1 < N£g < 
3. 

• Inverted mass hierarchy, mi ~ m 2 = m y 

In this model, m 3 is assumed to be the smallest 
mass, and m 3 < mi < m 2 . A peculiar feature of 
this model is that IV" r cannot be less than 2: when 
m 3 = 0 , the possible solutions allowed by the neu¬ 
trino oscillation experiments are m 2 — 0.055 eV 
and mi ~ 0.047 eV or m„ it0 t — 0.1 eV; thus, mi 
and m 2 must be very similar, giving TV" 1 ' — 2. 
Again, when m 3 is comparable to the mass dif¬ 
ference between m 2 and m 3 , all three masses are 
degenerate, TV"^ ~ 3. Hence, the “inverted mass 
hierarchy” model allows N£g to vary only in the 
limited parameter space, 2 < TV"| < 3. 

In the upper panel of Figure |3J we show APg for the 
normal mass hierarchy model (dashed line) and for the 
inverted mass hierarchy model (dotted line). One can see 
that the two models are indistinguishable (all masses are 
degenerate) for f v > 0 . 02 . 

How do we constrain 7V" r ? We measure m^tot from the 
overall suppression of power at small scales, as described 
in the previous section. Then we measure m v from 
the “break” of the power spectrum caused by the free- 
streaming scale, fcf s (z). The number of non-relativistic 
neutrinos is finally constrained as TV" r = m u> tot/m v , 
which tells us about the neutrino mass hierarchy. In the 
3rd to 5th columns in Table ITT1 we show the projected er¬ 
ror on TV" r from high-.z galaxy surveys combined with the 
Planck data, assuming the fiducial neutrino mass frac¬ 
tion of /„ = 0.05. As the error depends very much on 
the fiducial value of TV" r , we explore three different fidu¬ 
cial values, TV"g d = 1 , 2 and 3. As we have described 
above, however, the first two fiducial values, TV" dd = 1 
and 2 , are inconsistent with the neutrino oscillation ex¬ 
periments if /„ = 0.05 (all the masses must be nearly 
degenerate); thus, only the 5th column is actually real¬ 
istic if there is no sterile neutrino. We find that TV" r is 
going to be difficult to constrain: even when we com¬ 
bine all the high -z galaxy surveys, Gl, G2 and SG, the 
projected error is ATV" r = 0.6. This implies that it is 
not possible to discriminate between TV" r = 2 and 3 at 
more than 2 -a, while one can reject 7V" r = 1. When the 
fiducial value of /„ is small enough to allow for TV" r < 2, 
/„ = 0 . 01 , the constraints are too weak to be useful (see 
the 3rd to 5th columns in Table ITTT1) . 

As explained above, if the total neutrino mass is larger 
than 0.2 eV, Nffg ~ 3 is expected from the neutrino os¬ 
cillation experiments. The lower panel of Figure 0 shows 
the projected error on TV" 1 ' for the SG survey. A model 
with TV" r = 3 can be detected at more than ler level only 
if fu ~ 0.04 (m„, tot ~ 0.52 eV). Nevertheless, it should be 
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A full Planck (@ko = 0.05 Mpc 1 ) 
n s OLs 

@ k = 

^0,pivot (MpC ) 

&pivot 

n s 

a 3 

Planck alone 

0.0062 

0.0067 

- 

- 

- 

Gl 

0.0038 

0.0059 

0.030 

0.086 

0.035 

G2 

0.0037 

0.0057 

0.18 

0.018 

0.025 

SG 

0.0030 

0.0024 

0.48 

0.0033 

0.0070 


TABLE IV: Projected 68% errors on the parameters that characterize the shape of the primordial power spectrum, the tilt 
(n s ) and the running index (a s ). The left block is the same as the 6th and 7th columns in TablolTil which combines the galaxy 
survey data with the CMB data from Planck, while the right block shows the constraints when the CMB information on n s 
and a s are not used. The left block lists the constraints on n 3 and a 3 using ko = 0.05 Mpc -1 [see Eq. which was chosen 

such that the Planck data would yield the best constraints. The 1st column in the right block shows the pivot wavenumber at 
which the errors on n 3 and a s are uncorrelated for a given galaxy survey, and the 2nd and 3rd columns show the constraints 
on n 3 and a 3 at the pivot wavenumber, respectively. The space-based galaxy survey at 3.5 < z < 6.5, SG, on its own yields 
better constraints on n 3 and comparable constraints on a 3 compared to Planck alone, when evaluated at its pivot wavenumber. 
Note that fvfid = 0.05 and N^ d = 3 are assumed. 


noted that exploring the constraint on TV" r with future 
surveys is extremely important because any finding of 
a model with TV" r ^ 3 in this range of f u may provide 
valuable information for the existence of sterile neutrino 
or new physics. The general trend is that the error on 
_/V" r increases as m v decreases (i.e., the fiducial value 
of /„ decreases or 7V" r increases). This is because the 
neutrino free-streaming scale, fcf s , is proportional to m v : 
when nii, is too small, kf s will go out of the k range ac¬ 
cessible by galaxy surveys. A possible way to overcome 
this obstacle is to enlarge the survey volume which, in 
turn, can lower the minimum wavenumber sampled by 
the survey. The leftmost curve shows the projected error 
on TV" 1 ' for 7V"g d = 1 from a hypothetical full-sky SG 
survey at 3.5 < 2 < 6.5. We find that such a survey 
will be able to distinguish between two mass hierarchy 
models in principle. 


B. Shape of the Primordial Power Spectrum 

The amplitude of the primordial power spectrum ap¬ 
pears to be one of the most difficult parameters to mea¬ 
sure very accurately. Adding the galaxy survey does not 
help very much: the constraint on the amplitude im¬ 
proves only by a factor of 1.5 at most , even by combining 
all the data sets. (See the 9th column of Table El and 
mu) This is because the Planck experiment alone can 
provide sufficiently tight constraint on the amplitude: 
Planck allows us to break degeneracy between the am¬ 
plitude and the optical depth r by measuring the CMB 
polarization with high precision (the current accuracy of 
determining the amplitude, obtained from the WMAP, is 
about 10%). Adding galaxy surveys does not improve the 
accuracy of normalization due to the galaxy bias. The 
constraint on the amplitude could be further improved 
by adding the weak gravitational lensing data (e.g. see 
[531 tfiil [b 3|L which directly measures the dark matter 
distribution. Also, the lensing data are actually sensi¬ 



FIG. 5: The projected 68% error on n s (upper panel) and a s 
against fc max , as in Figure E 


tive to mass clustering in the non-linear regime, and thus 
would be complementary to the galaxy surveys probing 
the linear-regime fluctuations. 

The interesting parameter is the running spectral in¬ 
dex, a s . Actually G1 or G2 does not improve the con¬ 
straint on a s at all: the error shrinks merely by 20%; 
however, the space-based survey, such as CIP, provides 
a dramatic improvement over the Planck data, by a fac¬ 
tor of nearly 3. This indicates that SG alone is at least 
as powerful as Planck in terms of constraining a s . (We 
shall come back to this point below.) The driving force 
for this improvement is the value of the maximum us¬ 
able wavenumber for SG, fc max ~ 1 Mpc~\ which is sub¬ 
stantially greater than that for Gl, fc max ~ 0.2 Mpc -1 , 
and that for G2, fc max ~ 0.5 Mpc -1 . Our study there¬ 
fore indicates that one needs to push fc max at least up 
to fc max ~ 1 Mpc -1 in order to achieve a significant im¬ 
provement in the constraint on a s , for the survey parame¬ 
ters (the number density and bias parameters of sampled 
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FIG. 6: Left panel: Projected 68% error ellipses in the (n s , a s ) plane from Planck alone (the outermost contour), and the 
high -2 galaxy surveys combined with Planck. The dashed and dot-dashed contours are G1 (0.5 < z < 2) and G2 (2 < 2 < 4), 
respectively, while the dotted contour is SG (3.5 < 2 < 6.5). The highest redshift survey, SG, provides very tight constraints 
on the shape of the primordial power spectrum. The innermost shaded area shows the constraint from all the galaxy surveys 
and Planck combined. Right panel: Degradation in the constraints on n s and a s as a function of the non-relativistic neutrino 
contribution, /„, for SG plus Planck. The effect of non-relativistic neutrinos hardly affects the constraints on n s and a 3 . 
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FIG. 7: Projected 68% error ellipses in the {fv,n a ) (left panel) and (f„,a a ) (right panel) planes, respectively, for SG combined 
with Planck. Note that the fiducial values of f v ,Rd = 0.01 and AT“fi d = 1 are assumed. 




No CMB No Geometric Distortion No Redshift Space Distortion 


Survey 

U 

AC 

^m0 

n s 

a s 

U 

AC 

^mO 

Us 

OLs 

fu 

AC 

^mO 

n s 

a s 

G1 

0.19 

11 

0.037 

0.31 

0.22 

0.0052 

1.0 

0.0074 

0.0039 

0.0059 

0.017 

3.8 

0.0073 

0.0039 

0.0060 

G2 

0.082 

4.1 

0.023 

0.10 

0.060 

0.0043 

0.78 

0.0073 

0.0039 

0.0057 

0.015 

3.0 

0.0070 

0.0038 

0.0058 

SG 

0.044 

2.6 

0.0064 

0.057 

0.013 

0.0033 

0.64 0.0064 

0.0038 

0.0025 

0.011 

2.3 

0.0048 

0.0032 

0.0027 


TABLE V: Parameter degradation when some information is thrown away from the analysis. From left to right, CMB, the 
geometric distortion, or the redshift space distortion has been removed from the analysis. Note that = 0.05 and ACbd = 3 
are assumed. 
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V1N5 



V5N1 



V5N5 



fv 

jV nr 

n s 

a s 

fv 

N nr 

1 v V 

n s 

a s 

fv 

N? n a 

OL s 

Gl 

0.0038 

0.83(2.0) 

0.0037 

0.0059 

0.0033 

0.59(1.4) 

0.0036 

0.0057 

0.0029 

0.50(1.3) 0.0034 

0.0057 

G2 

0.0024 

0.71(1.8) 

0.0034 

0.0052 

0.0025 

0.59(1.4) 

0.0034 

0.0051 

0.0021 

0.46(1.2) 0.0028 

0.0040 

SG 

0.0018 

0.66(1.6) 

0.0026 

0.0018 

0.0017 

0.49(1.1) 

0.0023 

0.0018 

0.0016 

0.42(0.97) 0.0019 

0.0013 


TABLE VI: Projected 68% errors on the neutrino and power spectrum shape parameters for advanced survey parameters. 
“V1N5” has a factor of 5 larger number density of galaxies than the fiducial survey, “V5N1” has a factor of 5 larger survey 
volume, and “V5N5” has a factor of 5 increase in the both quantities. Note that f Vi fid = 0.05 and Nf T &d = 3 are assumed, while 
= 0.01 and A r " dd = 1 are assumed for Nf 1 in the parentheses. 


galaxies) we have considered. This can be also clearly 
found from Figure O 

On the other hand, the improvement on the tilt, n s , 
at fc 0 = 0.05 Mpc -1 is similar for Gl, G2, and SG: from 
a factor of 1.5 to 2. The interpretation of this result is, 
however, complicated by the fact that the actual con¬ 
straint depends very much on the value of the pivot 
scale, kg , at which n s is defined. The current value, 
ko = 0.05 Mpc -1 , was chosen such that the Planck data 
would provide the best constraint on n s . On the other 
hand, as the galaxy survey data probe fluctuations on 
the smaller spatial scales (larger fc), the optimal pivot 
wavenumber for the galaxy surveys should actually be 
larger than 0.05 Mpc -1 : the optimal pivot wavenum¬ 
bers for G2 and SG are fc’opivot = 0.18 and 0.48 Mpc -1 , 
respectively, where fco p i V ot was computed such that the 
covariance between n s (fcopivot) and a s should vanish at 
^Opivot and the two parameters would be statistically in¬ 
dependent (see [3l| for more discussion on this issue; see 
also [39f for the similar method for constraining the dark 
energy equation of state at pivot redshift). Table HVI lists 
^Opivot for Gl, G2 and SG, and the errors on n s (fcopivot) 
and a SiP ivot- Note that we do not use the CMB informa¬ 
tion on n s and a s to derive fcopivot (but include the CMB 
information on the other parameters). This Table there¬ 
fore basically shows how the galaxy survey data alone are 
sensitive to the shape of the primordial power spectrum. 
The striking one is SG: the errors on n s (fco P ivot) and a s 
are 0.0033 and 0.0070, respectively, which should be com¬ 
pared with those from the Planck data alone, 0.0062 and 
0.0067. Therefore, SG alone is at least as powerful as 
Planck, in terms of constraining the shape of the primor¬ 
dial power spectrum. 

The correct interpretation and summary of these re¬ 
sults is the following. The Planck data alone cannot con¬ 
strain the value of n s very well at small spatial scales 
(the uncertainty diverges at larger fco), and the galaxy 
survey data alone cannot do so at large spatial scales 
(the uncertainty diverges at smaller fco). However, when 
two data sets are combined, the accuracy in determin¬ 
ing n s becomes nearly uniform at all spatial scales, and 
the constraint becomes nearly independent of a particu¬ 
lar choice of fco . This is in fact a huge improvement of the 
situation, which may not be seen very clearly from just 
an improvement of n s defined at a particular fco- This 


is probably best represented by the constraint on a s we 
described above: the significant reduction in the uncer¬ 
tainty in a s for SG indicates that SG in combination with 
Planck has nearly uniform sensitivity to the shape of the 
primordial power spectrum from CMB to galaxy scales. 
This is exactly what one needs for improving constraints 
on inflationary models. 

The left panel of Figure [S] summarizes the constraints 
on a s and n s at fco = 0.05 Mpc -1 . This figure shows 
the projected error ellipses in the (n s , a s ) subspace. The 
overall improvement on the parameter constraint in the 
2D sub-space is quite impressive, and it clearly shows the 
importance of high -2 galaxy surveys for improving con¬ 
straints on the shape of the primordial power spectrum. 

Yet, one might be worried about the presence of non- 
relativistic neutrinos degrading the constraints, a s in 
particular, as the neutrinos might mimic the effect of 
a negative a s by suppressing the power more at smaller 
spatial scales. The right panel of Figurc[£j]basically shows 
that there is no need to worry: the constraints on n s and 
a s are hardly affected by the non-relativistic neutrinos 
for /„ = 0 — 0.05 (the errors on n s and a s are degraded 
less than by ~ 10%). This is because the inflation param¬ 
eters and the non-relativistic neutrinos change the shape 
of the galaxy power spectrum in modestly different ways, 
as explicitly demonstrated in the lower panel of Figure 
0 and the high -2 galaxy surveys are capable of discrimi¬ 
nating these effects. Figure [3 shows the error ellipses in 
the sub-spaces of (f v ,n s ) (left panel) and (f v ,a B ) (right 
panel). We find that the correlation between n s and f v is 
modest (the correlation coefficient defined by Eq. (fTTTll is 
0.55; see also Table lvnl) . while the correlation between a s 
and /„ is weaker. The most important result from this 
study is therefore that the 2D joint marginalized con¬ 
straint on inflationary parameters, n s and ct s , is hardly 
affected by the presence of non-relativistic neutrinos. 


C. Information from Geometric and Redshift 
Space Distortion 

In Table El we summarize what happens when we 
throw away some information from our analysis. Without 
CMB information, the errors on the neutrino parameters 
inflate significantly by more than an order of magnitude, 
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while the error on the matter density parameter, f2 m , is 
still comparable to or better than that from Planck alone. 
(Why this is so is explained in the next paragraph.) The 
errors on n s and a s at ko = 0.05 Mpc -1 also inflate. 

The geometric distortion effect helps to constrain 
from galaxy surveys alone. The radial distortion con¬ 
strains the expansion rate, H ( 2 ), while the transverse dis¬ 
tortion constrains the angular diameter distance, Da(z). 
Since H(z) and Da{z) have different dependences on Q ln 
and h for a flat universe, the distortion can determine O m 
and h simultaneously (I7 m and Q m h 2 for our parameter 
set). The galaxy surveys at higher redshifts benefit more 
from the geometric distortion effect. 

The redshift space distortion helps to constrain the 
neutrino parameters, by lifting degeneracy between the 
neutrino parameters and the galaxy bias. We demon¬ 
strate it in the third block of Table El The errors on 
the neutrino parameters increase up to nearly a factor of 
5 for SG, if we ignore the information from the redshift 
space distortion. The other parameters are not strongly 
affected. 


D. Variations with Survey Parameters 
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FIG. 8: The same as Figure^] but for advanced survey param¬ 
eters for SG combined with Planck, as in Table IyTI “V1N1” 
is the nominal survey parameters given in Table 0 Note that 
zero running, a a = 0, is assumed as a prior, and fu,Rd = 0.05 
and N™\ d = 3 are assumed. 


To guide the survey design, we show how parameter 
forecasts vary with two key survey parameters, the num¬ 
ber density of galaxies and the survey volume. Table EH 
shows how much one can reduce the projected errors on 
the key cosmological parameters by increasing the num¬ 
ber density of galaxies or the survey volume or both by a 
factor of 5. These advanced survey parameters are named 
as 

• V1N5: The survey volume is kept the same, while 
the number density of galaxies is increased by a 
factor of 5. 

• V5N1: The survey volume is increased by a factor 
of 5, while the number density of galaxies is kept 
the same. 

• V5N5: Both the survey volume and the number 
density of galaxies are increased by a factor of 5. 

Figure 0 shows how the error ellipses for f v and N™ will 
shrink for these advanced parameters. 

We find that the most effective way to improve deter¬ 
mination of the neutrino parameters, particularly iV" r , is 
to increase the survey volume. One may understand this 
from Figure^ — the information on the neutrino param¬ 
eters saturates at /c max ~ 1 Mpc -1 , and thus there is not 
much to gain by reducing the power spectrum errors at 
large k. 

On the other hand, one can still improve determina¬ 
tion of the parameters that determine the shape of the 
primordial power spectrum, n s and a s , by increasing the 
number density of galaxies. This is especially true for 
SG, which probes the largest k max . Therefore, one can 


achieve even better accuracies for constraining inflation¬ 
ary models by either increasing the survey volume or the 
number density of galaxies, in principle; however, in re¬ 
ality one is eventually going to be limited by our un¬ 
derstanding of the galaxy power spectrum in the weakly 
non-linear regime at large k. In other words, there is not 
much to gain by reducing the power spectrum errors at 
large fc, if the error is already as small as theoretical un¬ 
certainty in the modeling of galaxy power spectrum at 
the same k. Therefore, increasing the survey volume is 
probably the best way to improve accuracies of both the 
neutrino parameters and the inflationary parameters. 


VII. CONCLUSIONS AND DISCUSSIONS 

The non-relativistic neutrinos and the tilt and run¬ 
ning index of the primordial power spectrum cause scale- 
dependent modifications in the linear power spectrum 
probed by galaxy redshift surveys. We have shown that 
one can determine these parameters precisely by fully 
exploiting the two-dimensional information of galaxy 
clustering in angular and redshift directions from high- 
redshift galaxy surveys, when combined with the CMB 
data from the Planck experiment. The main results are 
summarized in Tables E] and Hill and may be graphically 
viewed in Figures 01 and 01 Our conclusions are two fold. 

The first conclusion is for the neutrino parameters. We 
have found that the future galaxy surveys with f l s = 300 
deg 2 can provide very tight constraints on the total neu¬ 
trino mass. The neutrino oscillation experiments have 
given the solid lower bound on the total neutrino mass, 
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m„.tot ~ 0.06 eV, which is actually larger than the pro¬ 
jected error on the total neutrino mass expected from 
the high-redshift galaxy surveys we have considered, up 
to by a factor of 2.5 for the space-based survey targeting 
galaxies at 3.5 < z < 6.5. If two neutrino species have 
nearly equal masses (the inverted hierarchy), then the 
lower bound from the neutrino oscillation experiments, 
m„,tot ~ 0.1 eV, is up to 4 times larger than the pro¬ 
jected errors of cosmological experiments. Overall, the 
high-redshift galaxy surveys combined with Planck al¬ 
low a positive detection of the total neutrino mass rather 
than the upper limit, improving the constraints on the 
total neutrino mass by a factor of 20-40 compared with 
the current cosmological constraints. The error on to i/) tot 
that we have found is smaller than that shown in the pre¬ 
vious work (e.g., see |hJ]) by a factor of 5 - 10, despite 
the fact that the survey volume we assumed is larger than 
theirs only by a factor 1.5 - 3. The main reason for the 
significant improvement is because our analysis exploits 
the full two-dimensional information in the galaxy power 
spectrum in redshift space. In particular, the redshift 
space distortion due to peculiar velocity field significantly 
helps improve the parameter determinations by breaking 
degeneracies between the cosmological parameters and 
the galaxy bias (also see Table [VJ. 

In addition, we have carefully investigated how one 
can use the future surveys to constrain the number of 
non-relativistic neutrino species, which should play an 
important role in resolving the neutrino mass hierarchy 
problem as well as the neutrino absolute mass scale, 
independently of the total neutrino mass. While we 
have found that the accuracy needed to discriminate be¬ 
tween two models with the same m [/j tot but different IV" r , 
cr(N™) £ 1, is going to be difficult to achieve for the 
nominal survey designs we considered, one may achieve 
the desired precision by enlarging the survey volume. It 
should be stressed here that it is extremely important 
to exploit independent constraints on m „ t tot and -/V" r 
from future cosmological observations, as any results un¬ 
expected from the point of view of the standard three- 
flavor neutrino model would imply anomalies in our un¬ 
derstanding of the neutrino physics and hints for the new 
physics. Needless to say, controlling the systematics in 
such observations is also a crucial issue. 

The second conclusion is for the shape of the primor¬ 
dial power spectrum. We have graphically summarized 
the expected performance of the future galaxy surveys 
for constraining the tilt, n s , and the running index, a s , 
of the primordial power spectrum in Figure 0 Com¬ 
pared with the constraints from the CMB data alone, 
the galaxy surveys we have considered can improve the 
determinations of n s and a s by a factor of 2 and 3, yield¬ 
ing a{n s ) ~ 0.003 and cr(a s ) ~ 0.002, respectively. The 
high-redshift galaxy surveys allow us to probe galaxy 
clustering in the linear regime down to smaller length 
scales than at low redshifts. It is also important to note 
that the galaxy survey and CMB are sensitive to the pri¬ 
mordial power spectrum shape at different k ranges, and 


therefore these two probes are nicely complementary to 
each other in terms of constraining the tilt and the run¬ 
ning index at different k. We have explicitly shown that 
the space-based galaxy survey (SG), such as CIP, has a 
similar level of precision for constraining n s and a s at 
the pivot scale fco p ivot ~ 0.5 Mpc -1 , to the Planck ex¬ 
periment at k 0 = 0.05 Mpc -1 (see Table ITvli . While it is 
sometimes argued that simple inflationary models should 
predict negligible amounts of running index which are out 
of reach of any experiments, our results show that the 
projected error on a s from the SG survey with five times 
more survey volume is actually as small as a s predicted 
by the simple model based on a massive, self-interacting 
scalar field, a s = —(0.8 — 1.2) x 10 -3 . Our results for n s 
and a s are still on the conservative side, as we have ig¬ 
nored the baryonic oscillations in the analysis. We have 
found that one can reduce the uncertainties in n s and 
a s by a factor of 2-3 by including the baryonic oscilla¬ 
tions — however, in that case one might also have to 
worry about the effect of curvature of the universe and 
dark energy properties. We report the results of the full 
analysis elsewhere. 

Throughout this paper, we have assumed the standard 
three -flavor neutrinos, some of which have become non- 
relativistic by the present epoch, as favored by neutrino 
oscillation experiments and the current theory of parti¬ 
cle physics. On the other hand, cosmological observa¬ 
tions such as CMB and the large-scale structure can also 
put independent constraints on the number of relativis¬ 
tic, weakly interacting particles just like neutrinos, as 
a change in the relativistic degrees of freedom directly 
affects the expansion rate of the universe during the ra¬ 
diation era. For example, an increase in the relativistic 
degrees of freedom delays the matter-radiation equality, 
to which the CMB and large-scale structure observables 
are sensitive. The number of relativistic particles (mi¬ 
nus photons) is conventionally expressed in terms of the 
effective number of neutrino species (i.e., the tempera¬ 
ture of additional relativistic species is assumed to be 
(4/11) 1 / 3 of the temperature of photons, just like or¬ 
dinary neutrinos). The current cosmological bound is 
N v = 4.2+]^ (95% C.L.) (68|, while the standard model 
predicts N v = 3.04 70] (also see (fitj for the constraint 
on the abundance and mass of the sterile neutrinos based 
on the recent cosmological data sets). Hence, it is inter¬ 
esting to explore how one can simultaneously constrain 
the number of relativistic degrees of freedom before the 
photon decoupling epoch as well as the number of non- 
relativistic particle species at low redshifts, _/V" r , using 
the future cosmological data sets. Properties of such 
weakly interacting particles are still difficult to measure 
exp er iment ally. 

It has been shown in the literature that the weak gravi¬ 
tational lensing miziEi or the number count of galaxy 
clusters j7J, IZa can also be a powerful probe of the cos¬ 
mological parameters. Different methods are sensitive to 
the structure formation at different ranges of redshifts 
and wavenumbers, and have different dependence on the 
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cosmological parameters. More importantly, they are 
subject to very different systematics. Hence, by combin¬ 
ing several methods including the galaxy surveys consid¬ 
ered in this paper, one can check for systematics inherent 
in one particular method. The combination of different 
methods may also reduce statistical errors on the cosmo¬ 
logical parameters. While it is worth exploring this issue 
carefully, we can easily imagine that, for the neutrino pa¬ 
rameters, the weak lensing or the cluster number count 
does not improve iV" r very much. To constrain 7V" r bet¬ 
ter one has to find a way to probe very large spatial scales, 
larger than the neutrino free-streaming scales; however, 
these methods probe only very small spatial scales where 
fluctuations have already become non-linear. (They may 
provide improvements in the determination of m„,tot-) 

There is a promising way to check for systematics and 
improve the parameter determinations using the galaxy 
survey alone. While we have been assuming that per¬ 
turbations are strictly linear and hence perfectly Gaus¬ 
sian, in reality small non-linearity always exists. The 
corrections to the power spectrum due to such small non¬ 
linearity may be calculated analytically usin g th e higher- 
order cosmological perturbation theory (see [76J for a re¬ 
view), which works extremely well at z > 2 (Komatsu 
and Jeong, in preparation). Using the same higher-order 
cosmological perturbation theory, one may also compute 
the higher-order statistics, such as the bispectrum, which 
is a very powerful tool to check for systematics due to 
non-linearity in matter clustering, redshift space distor¬ 
tion, and galaxy bias E30- Also, the bispectrum and 
power spectrum have different cosmological dependences 
|53L |73|; thus, it is naturally expected that combining 
the two would improve the determinations of cosmolog¬ 
ical parameters. The results of our investigation along 
these lines will be reported elsewhere. 

No matter how powerful the bispectrum could be in 
terms of checking for systematics, better models for non¬ 
linearity in redshift-space distortion and galaxy bias are 
definitely required for our projected errors to be actually 
realized. A good news is that we do not need a fully non¬ 
linear description of either component: we always restrict 
ourselves to the “weakly non-linear regime” where per¬ 
turbation theory should still be valid. Having an accurate 
model for the redshift space distortion in the weakly non¬ 
linear regime is important for the precise determination 
of the total neutrino mass, as the redshift distortion plays 
a major role in lifting the degeneracy between the galaxy 
bias and the matter power spectrum amplitude (see Ta¬ 
ble 0 what happens when the redshift space distortion is 
ignored). Recently, a sophisticated model of the distor¬ 
tion effect including the weak non-linear correction was 
developed in [7lJ based on the analytic method as well 
as the simulations. Likewise, it will be quite possible 
to develop a sufficiently accurate, well-calibrated model 
of the distortion effect at least on large length scales, 
based on adequate simulations. While we have employed 
a scale-independent linear bias throughout this paper, 
this model must break down even at weakly non-linear 


regime. Analytical as well as numerical [fill l62l l8ll| 
work has shown that deviations from a scale-independent 
bias do exist even on large spatial scales. This effect 
would also become particularly important for the preci¬ 
sion measurements of the baryonic oscillations for con¬ 
straining properties of dark energy. Therefore, careful 
and systematic investigations based on both analytical 
and numerical tools are needed to understand the real¬ 
istic effect of scale-dependent bias on estimates of the 
cosmological parameters. As we have mentioned already, 
information from the higher-order statistics in the galaxy 
clustering would be a powerful diagnosis tool to check for 
systematics due to non-linear bias. As perturbation the¬ 
ory predicts that the galaxy power spectrum and bispec¬ 
trum should depend on the galaxy bias differently, one 
can directly determine the galaxy bias and cosmological 
parameters simultaneously, by combining the two statis¬ 
tical quantities 0, 78| ,73, Ell Ei> ■ This method should 
also allow us to see a potential scale-dependent biasing 
effect from the epoch of reionization [f^] ■ 


Finally, let us comment on survey parameters. In order 
to make our discussion general, we have considered three 
hypothetical surveys which are different in their redshift 
coverage and the number density of targeted galaxies (see 
Sec. [3 for the survey definition). We have also explored 
how the parameter errors would change when the survey 
parameters are varied from the fiducial values (see Ta¬ 
ble [W] and Figure |HJ). Increasing a survey volume has a 
greater impact on the parameter errors compared with 
increasing the numbers of targeted galaxies for a given 
survey volume. We are hoping that our results provide 
useful information to help to define an optimal survey 
design to attain the desired accuracy on the parameter 
determinations, given the limited observational resources 
and budget. 
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APPENDIX A: PROPERTIES OF 
COSMOLOGICAL NEUTRINOS 


2. Non-relativistic Epoch 
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1. Mass Density 


The present-day mass density of non-relativistic neu¬ 
trinos is given by 


iv” 

Pu — ^ (Al) 

i=l 


where i runs over the neutrino species that are non- 
relativistic , and _/V" r is the number of non-relativistic 
neutrino species. We assume that some of the standard 
three active neutrinos are massive and thus 0 < N™ < 3. 
These neutrinos were in thermal equilibrium with other 
particles at early times until they decoupled from the 
primordial plasma slightly before electron-positron an¬ 
nihilation. Since they were still relativistic when they 
decoupled, their distribution function after decoupling is 
still given by that of a massless Fermion. After electron- 
positron annihilation, the temperature of photons be¬ 
came higher than the temperature of neutrinos by a fac¬ 
tor of (11/4) 1 / 3 . Thus, the neutrino number density of 
each species, n v , at a redshift relevant for a galaxy survey 
is given by the relativistic formula: 


n 


V 


3C(3) 3 

2 t r 2 " 


112(1 + z) 3 cm -3 , 


(A2) 


where £(3) ~ 1.202 and T v = (4/ll) 1 ' /3 T 70 (l + z), where 
T 7 o is the present-day photon temperature and we have 
assumed T 70 = 2.725 K [50f. Note that the number den¬ 
sity includes the contribution from anti-neutrinos as well, 
and does not depend on neutrino species, n v = The 
density parameter of non-relativistic neutrinos is thus 
given by 




8nGp™ 


8nGn 


IV” 


8H% 


i =1 


Ei 

94.1 h 2 eV' 


(A3) 


Since Vt v h 2 must be less than the density parameter of 
dark matter, fl c dmh 2 ~ 0.112, the total mass of non- 
relativistic neutrinos must satisfy the following cosmo¬ 
logical bound, 


Neutrinos became non-relativistic when the mean en¬ 
ergy per particle given by 

(A5) 

fell below the mass energy, m v y. The temperature at 
which a given neutrino species became non-relativistic, 
T”j, is thus given by 


rpr ir 
is,i 


180C(3) 

7tt 4 


m v ^i ~ 3680 



K. 


(A6) 


The redshift at which a given neutrino species became 
non-relativistic is 


1 + z nr;i ~ 1890 (^r) . (A7) 

As the current constraints from the galaxy power spec¬ 
trum at z ~ 0 already suggest < 1 eV, it is certain 
that neutrinos became non-relativistic during the matter- 
dominated era. The comoving wavenumber correspond¬ 
ing to the Hubble horizon size at z nr is thus given by 


, = H(z nI ,j) Q^hjl + Zn^j) 1 / 2 

" 1 + Znr.i 2998 M P C 

- 0.0145 (^) 1/2 O^/iMpc- 1 . (A8) 

Note that this value is smaller than that given in [3lj 
by a factor of ~ v/3.15. They assumed that neutrinos 
became non-relativistic when T u = m v , rather than T v = 
(180£(3)/77r 2 )m„ ~ mjz/3.15. 


3. Neutrino Free-streaming Scale 


Density perturbations of non-relativistic neutrinos 
grow only when the comoving wavenumber of perturba¬ 
tions is below the free-streaming scale, fcf s , given by 


^fs ,i(^) — 


H(z) 


2 (1 + z)a v ,i(z) ’ 


(A9) 


where a 2 
given in 


; Az ) is 
as 


is the velocity dispersion of neutrinos and 


d 3 p p 2 /m 2 


2 / \ _ / exp[p/T„(z)] + l 
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d 3 p 


15C(5 )T 2 (z) 

C(3) rn 2 i 


exp[p/T„(z)]+l 

15C(5) 2/3 T 2 (0)(1 + z) 2 

C(3) v 11 / 


, (A10) 


iv” 

i v< i < 10.5 eV. 

i=l 


(A4) 


where ((5) ~ 1.037. Hence, 


kf s ,i(z) 


0.677 

(1 + z) 1 / 2 



(All) 
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Neutrino density perturbations with k > k{ s ^ cannot 
grow because pressure gradient prevents neutrinos from 
collapsing gravitationally; thus, neutrinos are effectively 
smooth at k > fcfs.i; and the power spectrum of neu¬ 
trino perturbations is exponentially suppressed. Note 
that Eq. m is exactly the same as the Jeans scale in 
an expanding universe for collisional particles, if a v is 
replaced by the speed of sound. 


APPENDIX B: INFLATIONARY PREDICTIONS 


1. Generic Results 


Inflationary predictions are commonly expressed in 
terms of the shape of the primordial power spectrum of 
curvature perturbations in comoving gauge, 1Z: 

(ftk K') = (27r) 3 P K (k)S^(k - k'), (Bl) 

where 

1 , \ — l-t-n a -(-|a 3 ln(/c/fc 0 ) 

kj ■ (B2) 


k 3 Pn(k) 

2tt 2 


= 4 


Here, n s and a s are called the “tilt” and “running index” 
of the primordial power spectrum. These parameters are 
related to the shape of potential, V(fi), of an inflaton 
field, fi , a field which caused inflation, as Q 


n s - 1 


cx s 



2 M p 4 t 





(B4) 


where primes denote derivatives with respect to fi, and 
Mpi = ( 87 rG )" 1 / 2 = 2.4x 10 18 GeV is the reduced Planck 
mass. Successful inflationary models must yield suffi¬ 
ciently large number of e-folds for the expansion of the 
universe, N, before inflation ends at tend when fi rolls 
down on the potential to fi = fiend'- 

/•tend 1 r<t> Y 

N= dt H(t)~j7T d ^TFr ( B5 ) 

Jt M pl Jfand V 

which has to be at least as large as 50. This condition 
requires \n s — 1 | and |ct s | to be much less than unity, 
while exact values depend on specific inflationary models 
(i.e., the shape of V(fi)). Therefore, precision determina¬ 
tion of n s and a s is a very powerful tool for constraining 
inflationary models. 


2. Specific Examples 


An illustrative example of inflationary models is infla¬ 
tion caused by a massive, self-interacting real scalar field: 

V(fi) = ^mlfi 2 + fi A , (B 6 ) 

where is the mass of fi and A (> 0) is the coupling 
constant of self-interaction. The mass term and the in¬ 
teraction term equal when fi = fi c , where 


fic — 7T l(j) 



(B7) 


The mass term dominates when fi < fi c , while the inter¬ 
action term dominates when fi > fi c . One obtains 



Let us take the limit of mass-term driven inflation, 
fi <C fic- One finds 


N 


n s - 1 


a s 


fi 2 

4 M p y 

8M p 2 , _ 2 

fi 2 ~ N’ 

32M p 4 _ 2 

fi 4 ~ N 2 


(BH) 

(B12) 

(B13) 


For N = 50, n s = 0.96 and a s = —0.8 x 10 -3 . In the 
opposite limit, self-coupling driven inflation, fi fi c , one 
finds 


N 


n s 


1 


a s 


fi 2 

8M p V 

24M p 2 ! _ 3 

fi 2 ~ N’ 
192Af* _ 3 

fi 4 ~ N 2 


(B14) 

(B15) 

(B16) 


For N = 50, n s = 0.94 and a s = —1.2 x 10 -3 . These 
simple examples show that a precision of a(n s ) ~ 10^ 3 
is sufficient for discriminating between models, while 
cr(a s ) ~ 10 -3 may allow us to detect a s from simple 
(though not the simplest) inflationary models driven by 
a massive scalar field with self-coupling. 
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APPENDIX C: NORMALIZING POWER 
SPECTRUM 

While inflation predicts the power spectrum of 1Z, what 
we observe from galaxy redshift surveys is the power 
spectrum of matter density fluctuations, 5 m . In this sec¬ 
tion we derive the conversion from 1Z during inflation 
to 5 m at a particular redshift after the matter-radiation 
equality. 

Let us begin by writing the Poisson equation in Fourier 
space, 


which gives the power spectrum of 5 m , 

/ 2 k 2 \ 2 

P(k,a) = Pn(k)T\k)D 2 cbv (k,a), (C4) 

after the matter-radiation equality. The WMAP 
collaboration has determined the normalization of 
k 3 P-jz(k)/(2 tt 2 ) at k = k 0 = 0.05 Mpc -1 as 

= 2.95 x 10" 9 A, (C5) 


- k 2 ^ k (a) = 47rGp rn (a)5 m , k (a)a 2 = . 

a (Cl) 

where T is gravitational potential in the usual (Newto¬ 
nian) sense ( g 00 = — 1 + 24/). Cosmological perturbation 
theory relates 4/ after the matter-radiation equality to 7Z 
during inflation as 

T k («) = -\lZ^T(k) Dcb,j{k ' a \ (C2) 

where D cbl/ (k,a) is the linear growth factor of total 
matter perturbations including CDM, baryons and non- 
relativistic neutrinos, and T{k) is the linear transfer func¬ 
tion. Note that the transfer function and the growth 
rate are normalized such that T{k) —> 1 as k —> 0 and 
D cbu /a —> 1 as /c —> 0 during the matter era. (We as¬ 
sume that neutrinos became non-relativistic during the 
matter-dominated era.) Hence, 

2 k 2 

<5m,k(a) = „ 2n 7 Z k T(k)D cb „(k,a), (C3) 

O-tl Q ' im 


where A is a constant of order unity (Eq. [32] of Isfijb 
Putting these results together, we finally obtain the 
power spectrum of density perturbations normalized to 
WMAP: 


k 3 P(k, z) 
2-7T 2 


= 2.95 x 10~ y A 


2 k 2 


5H 2 n m 


xD 2 cb „(k,z)T 2 (k) 


-l+n,+|a s ln(fc/fc 0 ) 

(C6) 


APPENDIX D: CORRELATIONS BETWEEN 
PARAMETER ESTIMATES 


It is worth noting how the parameter estimations are 
correlated with each other for a given survey, which can 
be quantified by the correlation coefficients defined by 
Eq. 11911 . Table IVIII gives the correlation matrix for the 
parameters for SG combined with Planck. 
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